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We find a new correction to hyperfine splitting in the ground state of hydrogen atom in magnetic 
field. It can be interpreted as magnetic focusing of the wave function at the origin. The effect might 
be within the reach of experiment. 



m ■ 

o. 

cd ■ 



00 



Oh 



> 

o 
en 



The spectrum of hydrogen atom (HA) in strong niag- 
netic field (MF) was found long ago [1] and is presented 
in textbooks [2|. In recent years we are witnessing the 
rise of interest to this subject. This is probably due to the 
fact that huge MF up to eB - A^^-^ - 10^^ G has be- 
come a physical reality. Such field is created (for a short 
time) in heavy ion collisions at RHIC and LHC IS] . The 
field about four orders of magnitude less is anticipated 
to operate in magnetars |4J. Several interesting MF in- 
duced effects in QCD are under investigation now both 
from theoretical and experimantal sides |5|] . Among new 
results in physics of HA in MF necessary to mention the 
conclusion that in superstrong MF radiative corrections 
screen the Coulomb potential thus preventing the "fall to 
the center" phenomenon to occur [6|. In the present pa- 
per we discuss another MF induced effect, namely MF fo- 
cusing of hyperfine interaction. Experimentally it should 
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manifest itself via an additional shift of hyperfine states 
on top of the standard Zeeman splitting. 



We begin by introducing the units to be used and re- 
minding some basic equations. We put h — c — 1, 
a — e^ = 1/137, dimensionless MF is defined as iJ = 
B/Ba, Ba = m?e^ = 2.35 • 10^ G is the so-called atomic 
MF. At B — Ba the Bohr radius qb = {am)~^ becomes 
equal to the magnetic, or Landau, radius an = {eB)~^'^, 
the oscillator energy uj — eB/2m becomes equal to Ry- 
dberg energy Ry = mo? jl. In this system of units 
GeV^ = 1.45 • IQi^ G. 

The problem of HA in uniform MF is convenient to 
solve in cylindrical coordinates (p, z) using the London 
gauge A = iB X r, hence B is directed along the z-axis. 
The nonrelativistic Hamiltonian reads 
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where Aj^ is the transverse part of the Laplacian, [Ib = 
e/2m, e is the absolute value of the electron charge, cr = 
2se, Gz — ±1. The Schrodinger equation described by 
the Hamiltonian ([T|) does not allow the separation of the 
coordinates p and z. However in superstrong MF limit 
i? > 1 the "fast" MF variable p and the "slow" Coulomb 
variable z may be separated in the form of the adiabatic 
ansatz [l|,0]. 



Now we come to the subject of the paper. Hyperfine 
splitting (hfs) in the ground state of HA is measured to 
13 significant figures in frequency units [8|,|9| 



AEhfs = 1420.4057517667(9) MHz, 
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For 7? 3> 1 the dominant role is played by the lowest 
Landau level (LLL) with Up = 0, m — 0, —1, —2, ..., cr^ = 
— 1. For this state the energy of the oscillations in p- 
plane and the spin magnetic energy (IbB compensate 
each other. 



which amounts approximately to 5.9 • 10~^ eV. It corre- 
sponds to the 21cm line discovered in 1951 [10!| and since 
then thought to be primary tool in radioastronomy. The 
hfs can be found to lowest order in a from Breit magnetic 
dipole interaction 
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where gp = 2.79, ^in 
tion of Hhf gives 



^Ehf,^ — 



e/2mp. The first order perturba- 
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with n being the unit vector along the hne connecting dV 
and the origin where proton is placed. To proceed further 
we need the explicit expression for the wave function ^. 



There are three types of corrections to this expression: 
a)relativistic effects, b) QED, and c) nuclear structure. 
They have been thoroughly discussed in the literature - 
see e.g., p, 9]. 

With MF imposed equations (|3]) and (O experience 
important changes. In MF the HA takes the form of an 
elongated ellipsoid with r±_ ^ (_ff)~^/^, r^ ~ {\nH)~^ 
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Our calculations confirm this behavior, see Figl2] 
below. The corresponding expression for T-lhf can be de- 
rived making use of the Biot-Savart law [2|, |l2|. The 
operator Hhf has the form Hhf = —g^J'N{(Tp ■ B), where 
B is the MF created at the origin by the spin part of the 
electron current given by 



je = -A<bV*-^ X (Te 
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where function ^(p, z) is real and (^s- independent since 
we consider the ground state with l^ = 0. Next we have 
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n X j = -/is [V*2(n . ae) - cT,{n ■ V*^)] , (g) 



Attempts to find eigenvalues and eigenfunctions of 
Hamiltonian ([1]) have a long history - see [13| for a list of 
references. We use a variational method as many authors 
listed in [13| did. In certain features our approach bears 
a resemblance to that of IJ, |l5| . The wave function for 
the ground state is written as 
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where N = {■k'^ ' " r'^r ^ 



3/2j,2^j,^'j-i_ rpj^g ^^^ parameters r±_ and r^ 
are fitted at each value of H . The rationale for choosing 
'^Q in the form ([SJ is the following: a) it has a desired 
form of an elongated ellipsoid, b) it has an axial sym- 
metry and invariant under refiection with respect to p- 
plane, c) our calculations show that for H '^ \ the fitted 
wave function ([9]) is close to that obtained in [g, |l]| , and 
for < iJ < 1 the results are in agreement with very 
accurate calculations of several authors, e.g., [16|. 
The ground state energy is defined from 
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where Ho is obtained from ([T]) by removing the terms 
ujIz and iib^zB. Straightforward calculation leads to the 
following result for Eq 
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where 13 — r^/r^ < 1 for B > 0. Minimization of ([TT]) 
according to ([T0| yields r± and r^ as functions of H. For 
illustrative purposes consider two limiting cases: a)H — 
0, then uj = 0, r± ^ Tz, Eo = ima^/Sn ~ 0.85 Ry 
in line with [17|, b) free particle in MF, then we obtain 



(mw) 



Eo = LO. 



In Fig IT] we plot the energy Eo as a function of H 
in comparison with the results of [16J. The deviation 
from the elaborated calculation [16| does not exceed 15%. 
In Figl2] we display the radii r±^ and r^ as functions of 
H. This figure demonstrates how the deformation of the 
wave function with H proceeds. 



With the fitted wave function at our disposal we return 
to (O - dS]) and write the following expression for T-thf 
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The integrals can be evaluated analytically with the fol- 
lowing final result: 



Tihf - g^iBm [{FiiH) + F^(H)){(T, ■ (Tp)+ 



(13) 




FIG. 1. Plot of Eo (without spin contribution) vs. H. Solid 
curve - present calculation, dashed one from [I3] 




FIG. 2. The radii r± (solid hne) and r^ 
units as functions of H . 
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= r, and from (IT?| - p3|) one 
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As a result we return to Eqs. 
have 
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At i/ > 1 we 
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Fi - iJlniJ, F2 - \/^hi^ i/. 



(17) 



Equations p^ - dTB]) comprise the essence of the physi- 
cal process which can be called "Magnetic Focusing of 
Hyperfine Interaction" . MF compresses the HA thus in- 
creasing the wave function of the origin and giving rise 
to the tensor component. 
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The next task is to see how our results modify the 
standard Zeeman splitting effect. This question will be 
discussed at length in a forthcoming publication, while 
here we briefly outline some essential points. 

In MF the ground state of HA is splitted into four 
levels with their energies obtained by the diagonalization 
of the Hamiltonian 



ii'i^f = Uhf + l^B{(Te ■ B) - g^iN{o-p ■ B) 



(18) 



where Hhf is given by (J13p . Let us focus on the transi- 
tions between the states which at S = correspond to 
\a) = \S = 1,S^ = 0) and \b) == jS* = 0,5*^ = 0). From 
(fT8| one gets 



Ea-Eb^AEhfs\h^ + 
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where AEhjs is given by ^ and 7 — {Fi + F2)/2F. 
Without Magnetic Focusing 7 = 1 and the standard ex- 
pression is retrieved. The quantity of interest is the dif- 
ference Sv = V — vq with 1^0 corresponding to 7 — 1. 
Here we present semi-qualitative estimates of dv in the 
two limiting regimes of super-strong {H ^ 1) and weak 
{H <C 10^^) MF. Performing elementary calculations 
starting from (fT4|) . ([T5|) and (fT9|) we arrive at the fol- 
lowing results: 
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for iJ > 1 and 



for i?<a2^ 



di^ ~ AE, 



hfs 
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100 G. Evaluation of the quan- 
tity r^/r^ in the weak field limit requires accurate numer- 
ical calculations which will be presented in the forthcom- 
ing publication. We remind that the present hydrogen 
maser experiments are sensetive to the variations of the 
Zeeman splitting of the order of 1 mHz [18|. In Figl3]we 
show 5v in a rather wide interval of H. The growth of 
6v with H reflects the gradual deviation of HA from the 
spherical symmetry. 
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FIG. 3. The frequency shift Sv (see the text) vs. H. 



Magnetic Field Focusing considered here for the HA 
is a universal phenomenon important for any quantum 
system/reaction in presence of MF as soon as the wave 
function at the origin is an important parameter. In par- 
ticular, it leads to the modification of /3-decay rate in 
MF [19|. Another example is the spectrum of quark- 
antiquark system ;20!]. Interesting effects occur also in 
super-strong MF created in heavy-ion collisions. A few 
words are needed to add concerning related problems left 
beyond the scope of the present paper. Magnetic focusing 
in muonic hydrogen may be easier to observe experimen- 
tally [9]. For H ^ 1 another correction comes into play 
- proton can not be considered as infinitely heavy and 
problem becomes a two-body one [21[. 

The authors are grateful for many useful discussions 
and remarks to M.I. Vysotsky, S.I. Godunov, V.S. Popov 
and M.E. Eides. 
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